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Abstract. A leading candidate source of detectable gravitational waves is the inspiral and merger 
of pairs of stellar-mass compact objects. The advanced LIGO and advanced Virgo detectors will 
allow scientists to detect inspiral signals from more massive systems and at earlier times in the 
detector band, than with first generation detectors. The signal from a coalescence of two neutron 
stars is expected to stay in the sensitive band of advanced detectors for several minutes, thus allowing 
detection before the final coalescence of the system. 

In this work, the prospects of detecting inspiral signals prior to coalescence, and the possibility 
to derive a suitable sky area for source locations are investigated. As a large fraction of the signal is 
accumulated in the last ~ 10 seconds prior to coalescence, bandwidth and timing accuracy are largely 
accrued in the very last moments prior to coalescence. We use Monte Carlo techniques to estimate 
the accuracy of sky localization through networks of ground-based interferometers such as aLIGO 
and aVirgo. With the addition of the Japanese KAGRA detector, it is shown that the detection and 
triangulation before coalescence may be feasible. 



PACS numbers: 04.30.Tv, 04.80.Nn, 98.70.Rz 



1. Introduction 



In the near future advanced gravitational wave (GW) detectors, like advanced LIGO (aLIGO)[lJ, 
advanced Virgo (aVirgo)[2J and the Japanese KAGRAjf] [3] will be operative, allowing us to 
regularly observe gravitational wave signals. Nowadays there is a growing awareness that 
coincident observations of GW signals and electromagnetic (EM) counterparts will play a crucial 
role in the future of GW astronomy allowing us to distinguish between the cases of environmental 
rather than astrophysical origins of events in the data with a relatively low signal. Such 
observations will also be very important for astrophysics and cosmology. As an example, 
information regarding the supernova engine could be obtained by GW signals from core-collapse 
supernovae, GW signals coincident with gamma-ray bursts could help in determining the nature 
of the progenitor, and signals from a binary coalescence allow us to measure independently the 
distance and redshift to the source flU. This possibility will help us to determine the Hubble 
parameter to within a few percent EO. 

Because the major part of electromagnetic observatories are intrinsically directional, a possible 
approach to perform multi-messenger observations is to localize the position of the source from 
the gravitational wave observations as quickly and accurately as possible to use electromagnetic 
observatories to follow up the event. Then the first goal of GW astronomy would be to extract 
the sky location from the gravitational wave signal with enough accuracy. However, conversely 
to electromagnetic telescope, gravitational wave detectors are sensitive to most of the sky and it 
is not possible to obtain directional information of a short duration GW source from a single GW 
detector. Indeed the location of a GW signal is primarily obtained using triangulation methods 
based on the observed time delays of the signal at different detectors and, for this reason, a network 
of two or more detectors is crucial for multi-messenger astronomy. With the new era of network 
with three or more advanced detectors a better localization capability will be surely achieved. For 
example, a challenging project to follow up gravitational wave sources from initial GW detectors 
using electromagnetic telescopes is already underway B7l[8l. Localizing the GW sources with a 
good accuracy is not enough to perform multi-messenger observations, it must be done as rapidly 
as possible. Fortunately, the improvement of sensitivity in advanced detectors allow some GW 
signals from inspiral binaries to stay in the aLIGO sensitive band for several minutes, giving the 
possibility of an early detection. 

Regarding EM detection, a great number of electromagnetic observatories with a a field of 
view (FOV) of around ten square degree like Pan Starrs, Palomar Transient Factory, LOFAR, 
SkyMapper and many others (see J9l for an overview) are being specifically designed for transient 
phenomena. Obviously a good estimate of localization accuracy in the case of early detected GW 
signals will help us to understand the feasibility of the future follow-up searches. 
In this work we examine the ability of gravitational wave detectors to detect and localize transient 
signals before the coalescence. We perform a Monte Carlo simulation of sources distributed 
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uniformly in space for two different mass scenarios with the LIGO-Virgo (LV) and LIGO-Virgo- 
KAGRA (LVK) detector network. We compute the probabilities to detect a merger prior to 
coalescence, and the associated sky area. 

The layout of this paper is as follows. In section [2] the source population is described for this 
investigation, and section [3] describes how we calculate the SNR values of the simulated sources 
and how to obtain the sky area. Section [4] will summarize the results which are discussed in 
section |5l 



2. Source Population 

2.1. Waveforms 

In this work we will consider non spinning inspiral waveforms. Only the inspiral part of the 
waveform is used, which ends at the innermost stable circular orbit (ISCO) with frequency 
f isco = 6v /g^ GM • Here c is the speed of light, G is the gravitational constant and M is the total 
mass of the merger system. The frequency evolution of the waveform can be written as iflOl 

c 3 

where 

and rj = m ] ni2 - ]2 is the symmetric mass ratio. Eqs. (fB-pl) can be used to find the frequency of 



the system as a function of time. A measured signal h(t) can be written as 

h(t) = F+h+(jt) + F x h x (t) . (3) 

where h+(t) and h x (t) are the + and x polarization of the signal from the source, and F + and F x 
are the antenna response factors, which depend on the relative orientation between the location of 
the binary to the antenna ifTTl ITOl . As these factors depend on the location of the source relative to 
the detector, they are time dependent. Even for a low-mass binary system and a detector threshold 
of 10 Hz, the time the signal stays in the sensitive band of a detector is well below an hour. As 
most of the signal-to-noise ratio (SNR) is accumulated close to coalescence, the effects of Earth 
rotation can be neglected. 



2.2. Signal population 

We consider two populations of sources, one consisting of double neutron stars (NS-NS) with 
masses of 1.4M each, and one consisting of a neutron star with a black hole (NS-BH), with 
masses of 1.4M and 1OM , respectively. For the advanced detectors we assume two lower cutoff 
frequencies, one at / min = 30 Hz and a more optimistic at / min = 10 Hz. Systems consisting 
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of two black holes are not considered in this work, as the signal would coalesce at much smaller 
frequencies and the time spend in the sensitive detection band would be very short. Furthermore, 
it is unlikely that such systems would produce bright electromagnetic counterparts. 
The merger rate is assumed to follow the stellar birth rate in spiral galaxies, which can be 
determined by blue galaxy luminosities |[T2ll . However, there is reason to believe that the blue 
luminosity alone is not a good tracer of the merger rate, since a significant delay between the 
formation of the compact objects and their coalescence is expected. Therefore we assume a 
substantial fraction of merger to take place in old elliptical galaxies [fT3l [T4l . Since the range to 
which advanced detectors are sensitive is hundreds of Mpc, it is safe to assume a spatially uniform 
distribution of merger sources. 

3. Signal detection 

3.1. Single detector SNR calculation 

To find a signal in the data, we use standard matched filter algorithms lfi~5l [i~6ll . The SNR of the 
signal is 



where h*(f) is the complex conjugate of h(f) and Sh(f) is the power spectral density (PSD) 
function of the detector. The anticipated noise curves for advanced LIGO IfPTll and advanced 
Virgo lfi~8l are plotted in Fig. [T] As the Japanese KAmioka GRAvitational wave detector (KAGRA) 
might also be operational later in this decade, it is important to include this detector in this 
investigation. The power spectral density of KAGRA is available in text-form lfi~9l . which can 
be sufficiently fitted by the analytic expression 




(4) 



y/SM = 6.499 • 1(T 25 x 



(5) 



+ 9.72 • l(T 9 exp (-1.43 - 9.88x - 0.23x 2 ) 



+ 1.17exp (0.14 - 3.10x - 0.26x 2 ) 
+ 1.70 exp (0.14 + 1.09s - 0.013s 2 ) 

+ 1.25 exp (0.071 + 2.83x - 4.91x 2 ) 



where x = log(//100Hz). 

The advanced Virgo noise curve can be parametrized as 
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Figure 1. \J S(f) Power spectrum for Advanced LIGO, Advanced Virgo and KAGRA. 



y/SjJ) = 1-259 ■ 10- 



24 x 



(6) 



+ 0.07exp (-0.142 - 1.437s + 0.407a; 2 ) 
+ 3.10 exp (-0.466 - 1.043a; - 0.548a; 2 ) 
+ 0.40 exp (-0.304 + 2.896a; - 0.293a; 2 ) 

+ 0.09 exp (1.466 + 3.722a; - 0.984a; 2 ) 



where x = log(//300Hz). 

The central idea of this work is to estimate the accumulated SNR of a GW source before the 
coalescence of the system at the ISCO frequency. This is done by changing the upper limit on the 
integral in Eq. (01) 



p(/««) = a/4 l^^^fdf, (7) 



low 



and then using Eq. ([T]) to determine the maximum frequency / max as a function of time. The SNR 
can be expressed as ||20l 

M \ 5/6 /Mpc\ ^ /" /max /" ? / 3 

p = l.OO X 1U i - 

v 

where 



" = L56 x 10 " 19 [wj [Is )'"],. kuf' (8) 



^ g eo = y/F+0- - cos l2 ) 2 / 4 + F x cos L 2 (9) 
represents the response of an antenna in the direction of the source, i being the angle between the 
axis of orbital angular momentum of the binary and the line of sight to Earth. 
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Figure 2. Network antenna power pattern over the sky ^t.(-F+,fc + ^x,fc) f° r a three-detector 
network, consisting of the two aLIGO detectors and aVirgo. The network antenna power pattern 
describes the network relative sensitivity in different directions. Crosses represent 10000 detected 
sources 




Figure 3. Network antenna power pattern ^ fe +F X ^ ) over the sky for a four-detector network, 
consisting of the two aLIGO detectors, aVirgo and the Japanese KAGRA detector. The network 
antenna power pattern describes the network relative sensitivity in different directions. Crosses 
represent 10000 detected sources. 



3.2. Detector Network 

We now examine two different detector networks. The first network is made up of advanced 
LIGO and advanced Virgo (LV), while the second network includes the Japanese KAGRA detector 
(LVK). The expected PSD for advanced LIGO, advanced Virgo and KAGRA is used to calculate 
the estimated SNR of a source. 

A signal is considered to be detected by the network, if the recovered SNR exceeds a threshold of 
7 in each detector of the network ll2T1l . Although a reasonable detection of a signal can be made 
with two detectors alone, useful triangulation can only be done with three or more detector sites in 
the network. 
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Figure 4. SNR p and signal bandwidth 07 as a function of time prior to merger. Although a 
large fraction of the SNR is accumulated well before the coalescence, the bandwidth increases 
significantly at later times. The signal of a NS-NS binary enters the sensitive band at 52 seconds 
prior to coalescence (for a 30 Hz lower frequency cutoff), while it is 12 seconds for a NS-BH signal. 



The response of the three-detector network to sources in the sky is shown in Fig.[2| This plot shows 
the response of a given position above earth to the LV network. Each cross in this Figure represents 
a source in the sky that is being detected by the simulation, i.e. for which the recovered SNR is at 
least 7 in each detector. Fig. [3] shows the same result for the LVK network. The colorscale shows 
the network antenna power pattern ^2 k {F +> k + Fx,k) E2ll . where k goes over each detector. This 
is a generalization of the antenna power pattern for a network of detectors. It describes the relative 
sensitivity of the network in different directions. 



3.3. Localization of the source 

Given a list of detected sources, it is now possible to estimate the accuracy of localization for 
each source. For the purposes of this work, we assume the localization is achieved entirely 
by triangulation between the detectors although the individual amplitudes of the signal provide 
additional information at the ~10% level ||23ll . 

The accuracy of the observed time of coalescence a t depends on the SNR of the signal and the 
bandwidth Of of the detector's response. An approximation of this relationship can be found in 
Ref. El: 

°t = 7^— • (10) 

The timing accuracy is inversely proportional to the SNR p and the effective bandwidth a f of the 
source. 
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Figure 5. Number of detectable NS-NS signals (red) and NS-BH signals (green) as a function of 
time before coalescence. 



The approximations used to obtain Eq. ( 10) break down at low SNR, where second order effects 
become important. For a single detector SNR threshold of 7, this approximation differs from the 
exact calculation by ~10%, which is within the accuracy of this work (see right panel of Fig. 1 in 
Ref. 11231). 

3.3.1. Localization with three detectors There were several analytical Il24ll25l and numerical ||26l 
attempts to obtain a relation between time accuracy and the sky area. We will follow the analytic 
way laid out in Ref. [25] to obtain an approximate area given the observation of a binary signal. 
We approximate the area at a given confidence level CL as 

A(CL) ~ 2ir<j x a y [-]n(l-CL)]/cos9 , (11) 

where a x and a y denote the errors associated with the eigen-directions of the matrix describing the 
location accuracy (see Ref. [|24ll25ll for details), and 9 is the angle between the normal of the plane 
spanned by the three detectors and the direction to the source. 

Since the two LIGO detectors have comparable sensitivity and the light-travel time from each 
LIGO site to Virgo is about the same (27 ms), it is possible to simplify the analytical sky location 
expression of Ref. Il26ll by taking the timing accuracy 07 of the two LIGO detectors to be identical, 
but allowing the Virgo timing a v to differ. Using a median angle as cos 9 ~ 0.5, the 90% confidence 
level sky area A is 



e>7 



^^ 2odeg (o^sJ o^; ' • (i2) 



8 



Note, that with a network composed of three detectors, the sky area is degenerate since such a 
network cannot discriminate between a location above or below the plane spanned by the detectors. 
Since the error on the arrival time of a signal can be determined for each detector by using 
the bandwidth, it is important to note the properties of bandwidth before the coalescence. The 
dependence of the bandwidth and the SNR prior to coalescence is illustrated in Fig. |4j Although 
the SNR rises reasonably fast before the time of the merger, the bandwidth stays at very small 
values even very close to the merger because of the steep rise of the seismic wall in the power 
spectrum of the advanced detectors at lower frequencies (see Fig. [1}. 

3.3.2. Localization with four or more detectors The approximation described in the previous 
section to obtain the sky area in the case of a three detector network cannot be applied to networks 
with four or more detectors. Thus a different approach must be taken when dealing with a network 
with more detectors. In that approach, the whole sky is tiled into equidistant sky points, and the 
expected time delay between each detectors t p AB is calculated for each point P. The time residual 
for each detector combination is 



tAB - t^B 



(13) 



At AB 

where t A B is the observed light-travel time and At as is the error on the expected timing accuracy 
given as At AB = W a + a tB f° r detectors A and B. By computing the residual sum of squares 
over each detector combination combo, 

*T? SS = NJ2 **JZ*ab ; (14) 

, L±t>AB 
combo 

a probability can be assigned to each sky point P: 

p oc -e rss < . (15) 

After normalization of this quantities, all sky areas corresponding to the sky points are summed 
until the desired confidence level is reached. The resulting area A signal is the sky area where the 
source is located at the given confidence level. This method is described in more detail in Ref. [|23l . 



4. Results 



4.1. Single detector results 

Fig. [4] shows the evolution of the accumulated SNR and frequency bandwidth over the course of 
the inspiral for the NS-NS and the NS-BH case. In the case of a NS-BH binary, the SNR is 
accumulated much closer to coalescence compared to a NS-NS binary (see Table I). Since the 
gravitational wave signal from a NS-BH binary has smaller frequencies close to coalescence than 
a NS-NS binary, these systems reach the sensitive spectral part of a detector later, resulting in 
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Figure 6. Average sky area for a NS-NS system as a function of time. The error bars denote the 
standard deviation from the simulated sample. 



slower accumulation of SNR. Another interesting fact is the moment in which s signal enters 
the detector band. A signal from a NS-NS binary enters the sensitive band about 52 seconds 
prior coalescence, while it is only 12 seconds for a NS-BH system, given a frequency threshold 
of 30 Hz. If the frequency threshold is at 10 Hz, the NS-BH enters the detector band 215 
seconds before coalescence. Therefore, in order to detect more massive systems, a lower frequency 
threshold becomes crucial. Fig. [5] shows the fraction of detectable inspiral systems as time prior to 
coalescence, compared to the numbers detected at the time of coalescence. The fraction of NS-NS 
systems is much higher than that of NS-BH systems. At a time ~10 seconds prior to coalescence 
about 5% of the NS-NS systems are detectable, but only about 0.5% of the NS-BH systems. 
Early detection of a signal is only the first step to generate a notification for a telescope. The second 
step, the localization of the source, strongly depends on the timing accuracy of the signal. Fig. [4] 



shows that the bandwidth, and consequently the timing accuracy via Eq. ( 10), is largely accrued in 
the last moments prior to coalescence. In the next section we show that detecting and localizing a 
signal with a three-detector network alone is very challenging. Adding a fourth detector, like the 
Japanese KAGRA detector, results in dramatically increased detectability and location accuracy. 



Table 1. Fraction of the SNR as function of the time before coalescence for NS-NS and NS-BH 
systems, compared to the SNR accumulated at the time of coalescence. 
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4.2. Network Results 



In this section, the Monte Carlo simulation is applied to a detector network. A source is defined 
as detected, when the SNR in each detector i is p; > 7. This choice of threshold has been used in 
previous searches ll2"Tll . 

Future large-scale transient survey projects like the Large Synoptic Survey Telescope (LSST) 11271 . 
have an expected field-of-view (FOV) of ~10 square degrees. We estimate the probability p of 
observing the source of a GW signal inside that FOV by dividing the FOV by our estimated sky 
area. If the estimated sky area is smaller than 10 square degrees, the probability is set to 1: 



p — ) 6 /^signal "-^signal ^ -lw 



10deg 2 /A signa i iL4 signa i > lOdeg 
1 ifAignai < 10 deg 

For the three-detector LV network, the probabilities are listed in Table II and plotted in Fig. |7j The 
approximate average area for a signal detected prior to coalescence is shown in Fig. [6j Although 
the probability for catching the correct location at the time of coalescence is 90%, the probability 
decreases very fast and gets well below 1% for a signal detected <4 seconds prior to coalescence. 
Simulating the four-detector LVK network, the outcome is quite different. The additional detector 
in Japan introduces a network baseline, which is larger than for any detector pair of the LV network. 
Furthermore, the addition of a fourth detector eliminates the ambiguity of the sky area. However, 



because the approximation Eq. (12) can only be applied to a three detector network, we need to 



use the formalism for the general case as described in section 3.3.2 

When performing the Monte Carlo simulation with the LVK network, we found very encouraging 
results. Even for the most difficult case considered in this work, the detection and location of 
a NS-NS signal ~10 seconds prior to coalescence, the sky area is determined very precise. In 
98.5% of all cases, the sky area of the detected sources were localized to a sky area smaller 
than 10 square-degree. Since about 50 NS-NS sources a year are expected to be observed by 
advanced detectors lfT6l , about 2-3 sources may be detected and localized at ~10 seconds prior to 
coalescence. This opens very exciting prospects for science as the coalescence of a merger may be 
observed in coincident by GW detectors and optical telescopes ! In case of NS-BH systems those 
prospects are less encouraging, since the fraction that may be detected that early in a GW detector 
is only about 0.5%. However, the rate expectations for those systems are more uncertain, so it is 
very hard to make an adequate prediction. 

In order to get the most science of a coincident observation, an electromagnetic telescope must be 
pointed in the right direction in time. The whole data processing, data transfer and telescope 
pointing must not take more than 10 seconds. First, the data itself must be recorded at the 
detector sites, analyzed, and transferred to a central place, where the coincidence step is performed. 
This step might be accomplished within about 2 seconds [|28l . Second, if the analysis detects 
a coincident signal between the detectors, the sky area need to be computed and transferred to 
participating telescopes. This step may also be done within about 2 seconds. The third and last 
step, the actual telescope pointing, is most time consuming. For example, the pointing time of the 
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Figure 7. Probability of a GW-optical detection for a NS-NS system as a function of time before 
coalescence for the three detector network. 

Robotic Optical Transient Search Experiment (ROTSE) ||29l is about 4 seconds. To conclude, the 
total process from detecting a signal to the time a telescope is on target may be done in just under 
10 seconds. 

The sensitivity of a GW network is not uniform over the sky. The network is more sensitive to 
some sky regions than others, which can be seen in the sky maps shown in Fig.[2]and[3j The more 
sensitive regions on the sky are known a priori, and it is more likely a source if detected in this 
region. Telescopes could be pointed in such regions beforehand, or even an array of fixed mounted 
telescopes can be used to observe the interesting regions of the sky. Fig. [8] shows the fraction of 
the sky that need to be observed in order to make a detection of a signal at given probability. For 
example, if the GW network detects a signal from a GW source, a telescope array which covers 
40% of the whole sky should be able to catch this signal with a chance of ~80% . 

Table 2. Probabilities to have a 10-deg 2 telescope pointed to the correct location, in case the data 
analysis, the flow of information and the telescope's pointing time are negligible. 
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Figure 8. Fraction of detected sources as a function of the observed fraction of the sky. As the 
regions of the sky are known, in which the detector network can see farthest, a array of telescopes 
only need to cover some fraction of the sky. Covering 40% of the sky yield in a probability of around 
80% to have the source in the telescope's FOV. 

5. Discussion 

The work presented above shows that inspiral signals from the merger of two compact objects can 
be detected and localized before its coalescence. Although the majority of SNR is accumulated in 
the last moments before coalescence, an early detection and localization could become very likely 
with advanced detectors. Two GW detector networks have been considered, one consisting of the 
two aLIGO and the aVirgo detectors, and one network with the Japanese KAGRA in addition. 
The approximate SNR of simulated source is calculated using the expected spectral density noise 
curves. A source is considered detected, if the SNR is above 7 in each detector. The number of 
detectable sources decrease very steep at times prior to coalescence. Only ~3% of NS-NS are 
detectable at ~10 seconds prior to coalescence, while this number becomes ~0.6% for NS-BH 
sources. As the expected number of coalescences detectable with advanced LIGO/Virgo is around 
50 per year for NS-NS binaries, it might be possible to detect a few of such sources at about 10 
seconds before the coalescence. 

Having made a detection, the sky area must be computed from the observed signal. If the signal 
is observed with the 3-detector network, the sky areas are typically hundreds of square-degrees, 
which makes it very unlikely to have an array of telescopes pointed in time to the correct direction. 
However, when including KAGRA the prospects of early detection and location become much 
more likely. As KAGRA introduces a very long additional baseline to the network, avoids the 
ambiguity of a 3-detector network, and provides additional information to determine the sky area, 
it seems possible a single 10 square-degree telescope is enough to observe the sky area around the 
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coalescence. 

The question arises if there is enough time to process the GW data, and to send a notification 
to participating telescopes. With data-processing tools currently developed for advanced 
LIGO/Virgo, the data-analysis time can be reduced to virtually zero time. The time required to send 
the data between the detector sites and to the telescopes can also assumed to be sufficiently small. 
The most time consuming step seem to be the telescope pointing itself, but current telescopes are 
able to accomplish the pointing within a few seconds. Finally, optical observations during the 
coalescence of two compact objects seems possible. 

Even if a coincident GW/optical observation is possible only in a very few cases per year, this 
opens incredible scientific possibilities. It may become possible to directly test the connection of 
binary mergers with gamma-ray bursts, and to constrain models of their central engines. GRB 
outflow models can be tested, the geometry of the outflow, and tests of Lorentz invariance and 
tests of General Relativity in the strong field regime become possible. A complete overview of 
the implications of such a coincident observation is beyond the scope of this work, and subject of 
another investigation. 
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